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1.  INTRODUCTION 


This  progress  report  describes  the  work  performed  in  the  first  year  of  Project  2312/A5 
under  Grant  AFOSR-91-0122.  Die  aim  of  this  project  is  to  provide  a  theoretical  basis  for  the 
use  of  open  ended  coaxial  probes  for  dielectric  measurements.  More  specifically  we  had  three 
main  objectives: 

1.  To  model  the  fringing  fields  of  coaxial  sensors  in  lossy  media  in  order  to  determine 
the  extent  of  the  fringing  field  penetration  into  homogeneous  or  laminar  samples. 

2.  To  develop  an  admittance  model  for  the  probe  in  lossy  media  and  a  numerical 
technique  that  enable  its  use  for  dielectric  measurements. 

3.  To  evaluate  the  sources  of  error  associated  with  the  technique  and  to  develop  an  error 
model  to  estimate  the  sensitivity  of  the  measurements  at  each  frequency. 


Interest  in  the  use  of  open  ended  coaxial  probes  for  dielectric  measurements  stems  from 
its  suitability  to  a  variety  of  applications  in  material  science,  chemical,  biochemical  and 
biomedical  fields.  Its  acceptability  as  a  rigorous  measurement  technique  awaits  a  full 
understanding  of  the  theoretical  principles  involved.  We  are  currently  working  on  such  a  study 
and  will  be  presenting  our  preliminary  findings  in  this  report  The  theoretical  study  required  for 
the  development  of  an  admittance  model  has  now  been  completed  and  will  serve  as  a  basis  for 
the  study  of  field  penetration  in  lossy  dielectric.  The  assessment  of  the  errors  associated  with 
will  be  carried  out  in  the  very  near  future. 

Our  studies  show  that  the  simpler  admittance  models  more  commonly  used  for  dielectric 
measurements  can  be  derived  from  our  rigorous  formulation  provided  certain  simplifying 
assumptions  are  made.  The  full  solution  presented  in  this  report  provides  a  unifying  theory  for  a 
number  of  models  reported  in  the  literature  and  used  for  dielectric  measurements  with  various 

degrees  of  success^'^.We  have  studied  three  relatively  simpler  model  and  will  be  reporting  on 
their  range  of  applicability  and  their  limitations. 
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THEORETICAL  ANALYSIS 


2.1  Description  Of  The  Problem 


For  the  purpose  of  this  analysis  the  coaxial  probe  is  assumed  to  be  axisymmetric,  fitted 
with  a  ground  plane  that  is  large  by  comparison  to  its  diameter  and  that  in  addition  to  the 
principal  TEM  mode  only  TM^  evanescent  modes  that  preserve  cylindrical  symmetry  arc 
generated  at  the  interface. 

The  problem  is  therefore  best  described  in  the  cylindrical  coordinate  system  (p,  <J>,  z)  as 
illustrated  in  Figure  1  in  which  the  £j  refers  to  the  permittivity  of  the  dielectric  medium  inside 

the  line  and  outside  the  line ,  p  is  the  radial  distance  from  the  axis,  ({>  the  angular  displacement 

around  the  axis  and  z  the  displacement  along  the  axis.  In  the  first  instance  we  will  assume  that 
the  medium  outside  the  probe  is  uniform  and  occupies  the  entire  half-space  beyond  the  ground 
plate.  It  may  be  appropriate  to  mention  at  this  stage,  that  our  future  plans  include  making  a 
modification  to  the  formulation  to  treat  the  case  of  a  dielectric  material  of  finite  thickness,  backed 
by  an  infinite  perfectly  conducting  sheet  parallel  to  the  ground  plane. 


2.2  Expression  For  The  Electric  And  Magnetic  Fields  Inside  The  Line 

The  field  inside  the  line  are  a  superposition  of  the  forward  travelling  fundamental  TEM 
mode,  its  reflection  at  the  interface  and  the  TMq,,  evanescent  modes  generated  at  the 

discontinuity.  Their  derivation,  in  terms  of  a  potential  <t>,  from  Maxwell’s  equations  yields 
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EZ(P,Z)=X  f-^W^nCP) 

n=l  V  In  >/ 


(3) 


in  which 


r  is  the  reflection  coefficient  of  the  TEM  mode,  yn  is  the  propagation  constant  of  mode  n 

Yn  =  [  pS  -  k,2]10 


k!2  =  (02  Jl^i 

©  is  the  angular  frequency  and  rij  ,the  intrinsic  impedance  is 


for  n  >  0  pn  are  the  n  roots  of 

Jo(Pnb)  Y0(p„a)  -  J0(p„a)  Y0(pnb)  =  0 

where  J0,  Y0  arc  the  Bessel's  functions  of  the  first  and  second  kinds  of  order  zero,  a 
and  b  the  radii  of  the  coaxial  line. 


for  n=0  (TEM  mode)  pQ=  0 

For  each  n,  the  function  <Dn  is  defined  for  each  mode  as 


^n(p)  =  Jo(PnP)  Y0  (P„a)- J0  (p„a)  Y0  (pnp) 

with 

<Dn(a)  =  On(b)  =  0 

ao  and  an  are  amplitude  terms  normalised  to  the  dimensions  of  the  line  and  can  be 
expressed  in  terms  of  the  magnitude  of  the  electric  field  at  the  aperture  Ep*(p)  at  z=0  as 


and 


ao  =  |  E^(p)dp/(l+n 

ln(b/a) 

d<dpP)pEi(p)dp/[ 

dO„(p) 

.  dp  . 

2 

pdp 


with 
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2.3  Expression  For  The  Electric  And  Magnetic  Fields  outside  The  Line 


In  dielectric  medium  2  (z  >  0)  the  magnetic  field  is  expressed  in  terms  of  the  electric  field 
E-Kp)  at  the  aperture  as 


Hi(p)  =  i^l  10  Ep(P')  C0SW  ~  — p'dtp'dp' 


(4) 


where  p'  and  q>'  refer  to  the  coordinates  of  a  source  point  at  the  aperture  and  R  the 
distance  from  the  source  to  the  corresponding  field  point  which  at  z=0  gives 

R2  =  p2  +  p'2  -  2pp'  cos(<p'  -  <p) 

The  aperture  field  Ep*(p)  at  z=0  can  be  expressed  from  (2)  as 

E;<p)  =  ia0(l  +  r)  +  |;a„^^ 


(5) 


2,4  Expression  For  The  Terminating  Admittance 


Assuming  a  continuity  of  Ep*(p)  at  the  interface  the  continuity  of  which  across  the 
plane  z=0  allows  (1)  to  be  equated  to  H?-  (4)  which  gives 


^aod-n+i  E;(p')cos(<p'-9) ™ p'd<p'dp' 


(6) 


Substituting  for  and  introducing  the  terminating  admittance  Y*  of  a  coaxial  line  which 
is  defined  as 


where  Y0  is  the  characteristic  admittance  of  the  coaxial  line 

v  _  2* 

xo - 5~~ 

ln(— ) 
a 
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Equation  (6)  can  be  re-written  as 


.  f b  /*2Jt  -jk,R 

-i-fVp(p)dpY‘=y  &^(p)a+^2-  EpCpOcosCtp'-  <p)^— — p'dtp'dp' 

2rcpJ-  £  7„  dp  JaJ0  R 


Y‘=ly  f  f 2n Ep (p') 005(9' -  <p)  — p'dtp'dp' 

A£  y„  dp  2ji  JaJ0  R 


(7) 


where 


A  = 


1 

2ttp 


J*EJ(P)dp 


The  admittance  Y*  and  therefore  also  the  reflection  coefficient  F  can  be  calculated  from 
equation  7  using  a  appropriate  numerical  technique.  We  have  used  a  method  of  moment 
technique  (MOM)  for  a  full  solution  of  equation(7). 


3.  NUMERICAL  TECHNIQUES 

3,1  Full  Solution  Using  MOM 

The  MOM  solution  is  based  on  the  assumption  that  since  the  aperture  field  is  <p 
independent  then  within  a  small  concentric  ring  of  radius  pk  and  pk  ^  the  potential  Vk(p)  is 
constant  The  region  between  the  radii  of  the  coaxial  probe  a  and  b  is  subdivided  into  N 
intervals  with  boundaries  at  p0=a  at  which  Vk=l,  p^b.at  which  Vk=0. 

The  above  assumption  is  valid  if  dk  is  infinitesimally  small  which  require  N=°°.  A 
finite  number  of  rings  N  will  inevitably  introduce  a  truncation  error  To  overcome  the  problem 
we  have  introduced  a  quadratic  extrapolation  technique4  to  calculate  the  admittance  or  reflection 
coefficient  corresponding  to  N=«>  from  calculation  at  three  finite  values  of  N.  Figure  2  shows 
the  intermediate  and  extrapolated  results  of  reflection  coefficient  for  a  hypothetical  measurement 

configuration  of  a  coaxial  line  with  a=2. 333mm,  b=7.549mm  e1/£D=2.15  and  ej/E^lOO-jlOO  at 

1GHz.  This  rather  odd  combination  of  parameters  was  chosen  to  allow  for  comparison  with 
equivalent  values  in  reference  4.  The  agreement  was  found  to  be  of  the  order  of 
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10*6  and  was  deemed  satisfactory. 


3.2  Simplified  Variational  Methods 


The  equation  (7)  can  be  simplified  by  assuming  that  Ep*(p)  is  inversely  proportional  to 


p,  i.e. 


e;<p>=^ 


(8) 


where  Eq  is  a  constant  Incorporating  assumption  (8)  and  some  lengthy  but 
straightforward  algebraic  manipulations  (7)  becomes 


£0 

2ji 


jV 
J.  P 


Y‘=j<o e,X  „  [®nO>)-^ii(a)]J 


1 


Y„ 


rh- 


I  J^Eq 
2  n 


/-»* 


ra  /*2* 


a*' 


J 

aJ( 


a 

,-k,R 


d^n(p) 
dp 


pdp 


cos(p'-  <p) - dp'dp'dp 

o  R 


The  terms  of  the  n  summation  in  (9)  are  equal  to  zero  which  gives  an  aperture 
admittance  Y*  expressed  as 


(9) 


Y«  ^  mi 
[ln(b  /  a)]2 


j or- 


.-k,R 


cos(q>'-<p) - d<p'dp'dp 

R 


(10) 


This  expression  can  be  determined  by  expanding  the  exponential  term  into  an  infinite 
series  in  the  form 

Y*  =  coC0{e*C,  +  e^Qco2  +  e^/zQoP  +  e*3CJco4  +  e*7/2C6g>5  +  e^Qco6 

+  £*9/2 C.C07  +  £*50,0)8  +  e*i  i/2c,  0co?  +  ... )  (1 1) 

were  the  C  terms  are  constants  determined  by  the  dimensions  of  the  line  only 

We  have  used  expression  (1 1)  for  dielectric  measurements3  and  found  that,  in  practice, 
including  more  than  10  terms  had  no  effect  on  the  solution.  For  a  50  Q  Teflon-filled  coaxial 
probe  (3.5  mm  OD,  a  =  0.456  mm,  b  =  1.490  mm,  e^E,  =  2.1)  commonly  used  for  high 
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frequency  measurements,  equation  (11)  fails  to  converge  for  frequencies  in  excess  of  7  GHz  for 
lossy  dielectric. 

It  is  interesting  to  note  that  by  truncating  the  series  after  one  term  only  and  dropping  the 
higher  order  terms  of  ox  Equation  (11)  becomes5 


Y*=ii&(a+b)tE(m,-1, 


(12) 


where  E(m)  is  the  complete  elliptic  integral  of  the  second  kind  which  can  be 
approximated  by  a  polynomial,  and 


m  = 


4ab 

(a  +  b)2 


The  expression  (12)  can  thus  be  written  as 


Y*  =  jcoCe2 

and  is  equivalent  to  describing  the  admittance  of  the  probe  in  terms  of  a  frequency 
independent  lumped  capacitance  term.  Because  of  it  simplicity  formulations  equivalent  to 
equation  (12)  have  been  widely  used  for  dielectric  measurements. 


3,2  Alternative  Variational  Method 

The  simplifying  assumption  that  Ep*(p)  is  inversely  proportional  to  p,  is  implicit  in  (10) 
and  in  the  previous  section  we  have  shown  that  it  can  be  solved  by  expanding  the  exponential 
term.  An  alternative  solution  would  be  to  express  the  q>'  integral  in  terms  of  Bessel  functions 
which  gives*  after  changing  the  integration  parameters 


arg(^2-k|)V2  = 


0, 

7t 

~2' 


%>k2 

%<k2 
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Evaluating  the  p  integral  gives 


(13, 

The  integral  is  considered  separately  over  the  ranges  Oc£<k2  and  k2<4<°°.  The 
integrand  changes  from  imaginary  to  real  at  £=k2,  and  the  integration  over  these  two  ranges 
gives  the  aperture  conductance  G  and  susceptance  B  respectively.  For  0<£<k2,  the  variable  of 
integration  is  changed  to  £=k2  sin0  while  for  k2<^<«>,  the  present  infinite  integral  can  be 
expressed  as  a  finite  integral  of  integral  sines.  Given  that  for  the  coaxial  probe 


Y  =  G  +  jB 


The  solution  of  (13)  gives 


0.^ 


if; 


ln(tya)  J  sjn  q 


■  [j0(ko  Ve*bsin  0)  -  J0(ko  1  E*asin  "e)]2  d0 


(14a) 


B  = 


{  [2SiMeV  +  b2  -  2abcos  0))  -  2s{lK>^ e*a  sin(p) 
-  2Si|2ko^ e*b  sin  ^  d0 


(14b) 


where  k„  is  the  propagation  constant  in  free  space,  co/c,  e*  refers  to  medium  2 ;  J0  is  the 
Bessel  function  of  order  zero;  and  Si  is  the  sine  integral. 

The  admittance  model  described  in  (14a)  and  (14b)  is,  in  principle,  identical  to  that 
described  by  (1 1)  since  both  originate  from  (10).  The  only  difference  is  the  numerical  technique 
used  for  the  solution. 

A  summary  of  the  derivation  of  the  models  is  given  in  figure  3. 
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4.  PREDICTIONS  OF  ADMITTANCE  MODELS 


The  predictions  of  models  developed  in  the  previous  section  will  be  compared  with  each 
other  and  then  validated  against  measured  admittance  values  and  finally  used  to  perform 
dielectric  measurements  on  standard  solutions.  The  experimental  work  was  performed  on  an  HP 

8720  network  analyser  in  the  frequency  range  130  MHz  to  20  GHz  using  a  50  Q  Teflon-filled 

coaxial  probe  ( 3.5  mm  OD,  a  =  0.456  mm,  b  =  1.490  mm,  =  2.1,  21  mm  ground  plane). 

The  models  described  by  equations  10, 1 2  and  14  were  studied  in  details  in  a  previous 
report5.  Its  main  findings  were  as  follows: 

1.  For  the  probe  in  a  lossless  medium,  the  three  models  were  in  agreement  with  each 
other  and  with  the  measured  admittance  for  frequencies  up  to  20  GHz. 

2.  In  a  lossy  medium,  the  predictions  of  the  models  were  in  agreement  up  to  about  3 
GHz,  equations  10  and  14  were  in  agreement  for  frequencies  up  to  7  GHz. 

3.  Equations  14  gave  the  best  agreement  with  experimental  results  over  the  frequency 
range  studied. 

The  aim  of  this  report  is  to  describe  the  performance  of  the  newly  developed  model 
(equation  7)  and  to  compare  it  to  the  simpler  formulations.  Two  of  the  previous  models  (10  and 
12)  were  found  unsuitable  for  high  frequency  work  and  will  no  longer  be  mentioned.  The 
performance  of  the  new  model  will  be  compared  to  equations  14  only. 

Figures  4  to  6  show  a  comparison  between  the  predictions  of  equations  (7)  and  (14)  for 
the  real  and  imaginary  part  of  the  admittance  of  the  probe  in  standard  liquids  at  20°C  in  the 
frequency  range  100  MHz  to  20  GHz.  The  measurements  were  obtained  with  the  network 
analyser  calibrated  using  the  manufacturer’s  standards  prior  to  connecting  the  probe. 

The  predictions  of  both  models  exhibit  the  main  features  of  the  standard  liquids 
measured.  There  are,  however,  systematic  differences  betweem  the  two  models  and  between 
them  and  the  experimental  values.  It  important  to  mention  that  the  experimental  values  depends 
on  the  properties  of  the  standard  liquids  as  well  as  the  probe  and  connector.  To  illustrate  this 
point  measurements  on  water,  with  two  different  probes,  are  reported  in  Figure  7  showing 
clearly  that  the  experimental  values  fall  in  between  the  predictions  of  the  models. 
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5.  DIELECTRIC  MEASUREMENTS 


The  analysis  of  the  admittance  models  in  section  4  suggests  the  techniques  are  not 
suited  to  absolute  measurements  with  a  network  analyser  calibrated  using  factory  standards.  The 
models  assume  a  perfect  coaxial  probe  and  do  not  account  for  the  unwanted  reflections  of  the 
connector.  The  limitation  of  the  models  and  the  effect  of  the  probe  and  connector  can  be 
eliminated  if  either  model  is  used  as  part  of  a  calibration  at  the  end  of  the  probe.  This  is 
illustrated  in  Figure  8  were  the  predictions  of  equations  14  a  &b  are  compared  to  the  measured 
values. 


Using  the  end  of  probe  calibration  technique  both  models  were  successfully  used  for  the 
measurement  of  the  dielectric  properties  of  standard  liquids.  Figure  9  shows  the  agreement 
between  the  measurements  using  the  two  models  and  the  literature  values  for  formamide  at 
20°C. 


Experimentally,  the  main  difference  between  the  models  is  in  the  computation  time. 
Using  an  HP345  work  station,  it  takes  less  than  a  minute  to  compute  the  dielectric  properties  at 
21  frequency  points  using  equations  14.  It  takes  approximately  3  minutes  per  frequency  point 
using  the  model  described  by  equation  7. 


6  SUMMARY  AND  FUTURE  WORK 
The  aims  of  this  project  are: 

1.  Development  of  admittance  models  incorporating  the  excitation  of  TEM  and 
higher  order  modes  at  the  probe-sample  interface. 

2.  Coaxial  probe  field  penetration  in  lossy  media; 

3.  Error  analysis  and  model  evaluation 


Most  of  the  work  related  to  1  and  2  above  have  now  been  completed.  We  plan  to  cany 
out  the  error  analysis  in  the  near  future  before  embarking  on  a  series  of  dielectric  measurements 
using  the  new  techniques.  The  material  under  consideration  are  liquid  crystals  and  conductive 
plastics. 

Finally,  we  plan  to  write  two  or  more  publications  based  on  the  material  produced 
during  the  project 
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Amplitude 


Figure  1. Cylindrical  coordinates  of  a  coaxial  probe  probe  on  inner 
and  outer  radii  a  and  b  respectively. 
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Figure  2  Illustration  of  the  extrapolation  technique 
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Figure  3:  Derivation  of  admittance  models  for  a  coaxial  probe 


Conductance  of  probe  In  water 


— 


Equation  7 
Equation  14b 
Experimental 


Frequency  /Hz 


Susceptance  of  the  probe  In  water 


•  ■jAI 


“““  Equation  7 
Equation  14b 
B  Experimental 


6 


I»(Y) 


Fr«qu#ncy  /  CHx 

Figures  7 :  (a)  Real  and  (b)  imaginary  parts  of  die  admittance  of  8  coaxial  probe  in  water  at 
23.8°G  The  joined  symbols  are  experimental  measurements  and  the  solid  lines  die  predictions 

of  the  two  models 
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Figure  8  Measured  and  calculated  values  of  the  complex  admittance  of  the  probe  in  contact  with 
water  at  20  °C  and  in  the  frequency  range  100  MHz  to  20  GHz 

PERMITTIVITY  OF  FORMAMIDE  AT  20  C 


Figure  9:  Measured  (symbols)  and  literature  values  of  the  complex  permittivity  of  formamide  at 
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